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Abstract 

In this paper, results of 2. 5- dimensional magnetohydrodynamical simulations are 
reported for the magnetic reconnection of non-perfectly antiparallel magnetic fields. 
The magnetic field has a component perpendicular to the computational plane, that 
is, guide field. The angle 6 between magnetic field lines in two half regions is a 
key parameter in our simulations whereas the initial distribution of the plasma is 
assumed to be simple; density and pressure are uniform except for the current sheet 
region. Alfven waves are generated at the reconnection point and propagate along 
the reconnected field line. The energy fluxes of the Alfven waves and magneto- 
acoustic waves (slow mode and fast mode) generated by the magnetic reconnection 
are measured. Each flux shows the similar time evolution independent of 6. The 
percentage of the energies (time integral of energy fluxes) carried by the Alfven waves 
and magneto-acoustic waves to the released magnetic energy are calculated. The 
Alfven waves carry 38.9%, 36.0%, and 29.5% of the released magnetic energy at the 
maximum {6 = 80°) in the case of /3 = 0.1, 1, and 20 respectively, where /3 is the 
plasma /3 (the ratio of gas pressure to magnetic pressure). The magneto-acoustic 
waves carry 16.2% {9 = 70°), 25.9% {9 = 60°), and 75.0% {9 = 180°) of the energy at 
the maximum. Implications of these results for solar coronal heating and acceleration 
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of high-speed solar wind are discussed. 
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1. Introduction 

The heating mechanism of the solar corona is one of the most mysterious issue in astro- 
physics (Aschwanden 2004). The structure of the solar atmosphere consists of the photosphere, 
the temperature of which is about or more than six thousand degrees, the chromosphere, the 
temperature of which is up to a few ten thousand degrees, the transition region, in which the 
temperature increases rapidly, and the corona, the temperature of which reaches a few mil- 
lion degrees. To maintain such a high temperature in the corona in spite of cooling by heat 
conduction and radiative losses, a continuous supply of thermal energy is necessary. 

X-ray observations from early space experiments (e.g., Skylab) have shown that the 
corona is not uniform and consists of many bright loops. Poletto et al. (1975) demonstrated a 
correspondence between enhanced X-ray emission and a magnetic loop. It is therefore suggested 
that magnetic activity is related to the heating of the solar corona. Related to the magnetic 
activity, there are two promising models for coronal heating (see reviews by Klimchuk 2006, 
Erdelyi & Ballai 2007 and also references therein). One is heating by the dissipation of Alfven 
waves ^ that propagate in the magnetic flux tubes (Alfven 1947; HoUweg 1973, 1981, 1984, 1986; 
Uchida & Kaburaki 1974; Wentzel 1974; McKenzie, Banaszkiewicz, & Axford 1995; Axford et 
al. 1999). There are many proposed mechanisms for the dissipation of Alfven waves, such as 
mode conversion (see the next paragraph), resonant absorption (e.g., lonson 1978; Poedts et al. 
1989; Erdelyi & Goossens 1995; Ruderman et al. 1997), phase mixing (e.g, Heyvaerts & Priest 
1983), or magnetohydrodynamic (MHD) turbulence (e.g., Inverarity & Priest 1995; Matthaeus 
et al. 1999). The other is heating by many small-scale flares, i.e., nanoflares triggered by 
magnetic reconnection (Parker 1988). Observationally the occurrence frequency of microflares 
and nanoflares, A^, has been found to be dN/ dW oc VT"", where W is the total flare energy. The 
power-law index, a, ranges from 1.5 to 1.8 (Hudson 1991; Shimizu 1995; Shimojo & Shibata 
1999; Aschwanden & Parnell 2002) and nanoflares can not account for coronal heating if a is 
less than 2. However, a value of a larger than 2 has also been found (Krucker & Benz 1998; 
Parnell & Jupp 2000). The deflnitive conclusion has not yet been obtained from observations. 

As an origin of Alfven waves, Kudoh & Shibata (1999) considered a photospheric random 
motion propagating along an open magnetic flux tube in the solar atomosphere, and performed 
1.5-dimensional (1.5D, i.e., torsional motion is allowed) MHD simulations for solar spicule 
formation and the heating of the corona. It was shown that Alfven waves transport sufficient 

^ Other types of waves (acoustic, slow-mode and fast-mode waves) are strongly damped or reflected at the 
steep density and temperature gradients of chromosphere and transition region. 
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energy flux into the corona to account for its lieating, by extending tlie work by Hollweg, 
Jackson, & Galloway (1982) (see also Saito, Kudoh, & Shibata 2001). Moriyasu et al. (2004) 
performed 1.5D MHD simulations of the propagation of nonlinear Alfven waves along a closed 
magnetic loop including heat conduction and radiative cooling. They found that the corona is 
episodically heated by fast- and slow-mode MHD shocks generated by nonlinear Alfven waves 
via nonlinear mode-coupling. It was also found that the time variation of the simulated extreme- 
ultraviolet (EUV) and X-ray intensities is quite similar to the observed one. They concluded 
that the observed nanofiares may not be a result of reconnection but may be due to nonlinear 
Alfven waves. Subsequently, Antolin et al. (2008) discussed the observational signatures of 
the power-law indexes and coronal heating mechanisms (Alfven waves and nanofiares) by using 
1.5D MHD model of Moriyasu et al. (2004). They found that Alfven heating and nanofiare 
heating exhibit different power-law indexes (see also Antolin & Shibata 2010). 

However, since Alfven waves can be generated by magnetic reconnection unless the 
reconnection takes place in perfectly antiparallel magnetic fields, the nanofiare model may not 
be very different from the Alfven wave model. Yokoyama & Shibata (1995) modelled jets (X- 
ray or EUV jets and serges observed with Ha in the chromosphere) by performing a resistive 
two-dimensional MHD simulation of the magnetic reconnection occurring in the current sheet 
between emerging magnetic flux and overlying pre-existing coronal magnetic fields. Recent 
Hinode (Kosugi et al. 2007) observation revealed that jets are ubiquitous in the chromosphere 
(Shibata et al. 2007). Nishizuka et al. (2008) proved that the jets on the basis of this model are 
quantitatively corresponding to the multiwavelength jets (see also Cirtain et al. 2007 and Liu 
et al. 2009) observed with Hmode and TRACE (Handy et al. 1999). Yokoyama (1998) found 
that the ratio of energy of Alfven waves to the energy released by the magnetic reconnection 
(the model of Yokoyama & Shibata 1995) is nearly equal to 3%. In this model, there is a 
shear between emerging magnetic fields and coronal magnetic fields, kinks are produced by the 
magnetic reconnection and propagate away as Alfven waves. 

Takeuchi & Shibata (2001) performed 2.5-dimensional (2.5D) MHD simulations of the 
photospheric magnetic reconnection caused by convection, and found that the energy flux 
of Alfven waves ^ is enough to explain both coronal heating and spicule production. The 
generation of Alfven waves through magnetic reconnection has been also discussed by Parker 
(1991), Axford et al. (1999), and Sturrock (1999) in the context of coronal hole heating and 
solar wind acceleration. 

The upward acoustic waves in the flux tube are expected to be transported as other 
waves (e.g., slow- mode waves) after the formation of the shocks. Slow- mode waves are also 
thought to contribute to coronal heating though the contribution is only to the relatively low 
corona because of their compressibility. Takeuchi & Shibata (2001) measured the energy flux 
of slow-mode waves and found that it was ten times larger than that of Alfven waves. Suzuki 

^ The perpendicular magnetic field injected by magnetic reconnection propagates as Alfven waves. 
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(2002) discussed the possibility of coronal heating and the acceleration of the low-speed solar 
wind by slow-mode waves. Suzuki (2004) developed his study by including fast-mode waves 
(linearly polarized Alfven waves) and showed that slow-mode waves contribute to the heating of 
the low corona and the acceleration of the low-speed solar wind, while linearly polarized Alfven 
waves contribute to the heating of the outer corona and the acceleration of the high-speed solar 
wind. 

As the mention above, waves are created by torsional motions as recently observed (e.g., 
Bonet et al. 2008) in the photosphere or magnetic reconnections. They propagate into the 
corona and disipate their energy through linear and non-linear mechanisms. In particular, 
Alfven waves can be generated by such drivers, for the first time, by Jess et al. (2009). From 
Hinode data, De Pontieu et al. (2007) estimated the energy flux carried by transversal oscilla- 
tions generated by spicules and compared with radiative MHD simulations by more realistically 
extending Kudoh & Shibata (1999). They indicated that the calculated energy flux is enough 
to heat the quiet corona and to accelerate the high-speed solar wind. Okamoto et al. (2007) 
also estimated the energy flux to be 2 x 10^ ergs cm~^ s~^ propagating on coronal magnetic 
flelds. These reports (see also Tomczyk et al. 2007), however, have considerable argument 
about what is Alfven waves. Erdelyi and Fedun (2007) and Van Doorsselaere et al. (2008) 
argued that these oscillations were likely to be kink oscillations from observed behavior (see 
also Goossens et al. 2009; Taroyan & Erdelyi 2009). 

It is thought that the magnetic reconnection causes the solar flare. Particle acceleration 
takes place associated with solar flare. The mechanism of the particle acceleration, however, 
is not made clear. Alfven waves could contribute to the acceleration of ions through cyclotron 
resonance (see, e.g.. Miller 2000, and references therein). The generation of Alfven waves by 
the magnetic reconnection is a very interesting research topic also from such a point of view. 

In this paper, we present the results of 2.5D MHD simulations of the Alfven wave gen- 
eration by the reconnection of non-perfectly antiparallel magnetic flelds. The initial magnetic 
flelds have a shear, i.e., the magnetic fleld has a component perpendicular to the computational 
plane. This component can not contribute to the reconnection and is so-called guide fleld. The 
magnetic reconnection in this geometry was analytically studied by Petschek & Thorne (1967). 
The angle between the magnetic fleld lines in two half regions is a parameter, 6, while the 
initial distribution of the plasma is assumed to be simple. The energy fluxes of Alfven waves 
and magneto-acoustic waves are measured. In section 2 we describe the numerical method and 
model. In section 3 we show the results of the simulations. Finally a discussion and summary 
are given in section 4. 
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2. Numerical Simulations 



2.1. Assumptions and Basic Equations 

In this paper, we assume that the distributions of the physical quantities are not depen- 
dent on the ^-coordinate but that the vector fields have the 2;-component (i.e., 2.5D approxi- 
mation). We solve the following resistive MHD equations numerically: 

^ + {v.V)p = -p{V-v), (1) 

dv 11 

+ . V)-*; = —Vp+ — U X B) , (2) 
at p Anp 

^ + . V)|> = -ip{V ■ ^) + (7 - , (3) 

^^-VxE. (4) 

E = r]j — V X B, (5) 

and 

j = VxB, (6) 

where p,p,v are the density, pressure, and velocity of the gas, B,E,j are the magnetic field, 
electric field, and current density, and 7, 77 represent the ratio of specific heats and electric 
resistivity. 

2.2. Initial Conditions and Boundary Conditions 

We consider the situation that the physical quantities are uniform far from the current 



sheet. Their typical quantities are po = 1; Po = 1/7; and Bq = \Bo\ = y Sir / {13'y) , where Bq is 
the initial magnetic field. /3 = Supq/Bq^ is the plasma /3 (the ratio of gas pressure to magnetic 
pressure). Initially the magnetic pressure gradient force balances with the gas pressure gradient 



force. The sound speed, Cs = y'yp/p, is unity in the whole domain of the simulation box. 
The initial distributions of the physical quantities are as follows: 
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P = 7P, (11) 



where 




(12) 



Lq is the thickness of the initial current sheet. In this study, we have adopted 7 = 5/3 and 



Figure 1 shows the coordinate system [xy-plane and yz-plane), the initial magnetic field 
configuration, and the definition of 6 and (the directions of the wave vector and B± discussed 
in section 3.1 are also displayed). When 6 is equal to 180° the magnetic fields are perfectly 
antiparallel on both sides of a line where a: = 0. We investigated the cases where 6 is from 
180° to 10°. The value of the plasma /3 is also a variable parameter and equal to 0.1, 1, or 20. 
All physical quantities are normalized by their typical values. This means that the velocity is 



normalized by the initial sound speed, Cso = ylPo/po- Time is normalized by to = Lq/Cso- 



The number of grid points in the simulations is 603 x 1003. The grid spacing is uniform 
within — 5Lo < a; < 5Lq and — 20Lo < y < 20Lo. The constant grid spacing in the a;-direction 
is equal to 0.025Lo and that in the y-direction is O.OSLq. The grid spacing in the non-uniform 
region is slowly stretched by an increment of 2 % at each grid step (e.g., |Axj+i| = 1.02|Axj|). 
The size of the computational domain is — 13.1Lo < x < 13.1Lo, — 36.3Lo < y < 36.3Lo. 

At all the boundaries of the simulation box the periodic boundary condition is assumed. 
The total energy is therefore conserved. It can be said that the decrease of the magnetic energy 
in the simulation box is equal to the energy released by the magnetic reconnection. 

2.3. Resistivity Model and Numerical Method 

In order to initiate the magnetic reconnection, we assume a localized resistivity near the 
origin, {x,y) = (0,0), as 



where r = (x^ + y^) , = Lq, and rjo = O.ICsqLq. The magnetic Reynolds number is -Rm = 
VAoLy/rjQ, where v^o is the initial Alfven velocity outside the current sheet and Ly{= i/max ~ 
l/min = 72.6Lo) is the size of the computational domain in the y-direction. The values of the 
plasma /3 and i?m in our simulations are summarized in Table 1. We here mention that this 
non-dimensional parameter, R^, is much smaller than that of typical solar corona plasma 
(i?m ~ 10^^). However the numerical treatment of -Rm ~ 10^^ is very difficult in current-day 
computer resources. 

The numerical computations were carried out by the Rational CIP (Cubic interpolated 
profile) method (Yabe & Aoki 1991; Xiao, Yabe & Ito 1996) combined with the MOC-CT 
method (Evans & Hawley 1988; Stone & Norman 1992). The magnetic induction equation 
was solved by the MOC-CT and the other equations were solved by the CIP (e.g., Kudoh, 
Matsumoto & Shibata 1998; Kigure & Shibata 2005; Takahashi et al. 2009). 



'0 — 



1. 





(13) 
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3. Results of Numerical Simulations 



3.1. Time Evolution 

First we describe the time evolution of the system only in the case of plasma /3 = 0.1. 
Figures 2a and 2b show the distribution of the logarithmic pressure in the case of 6' = 180°. 
The magnetic fields reconnect in the diffusion region, which is around {x,y) = (0,0) and the 
magnetic energy is converted to thermal energy. The plasma is accelerated in the ^/-direction by 
the magnetic tension of the reconnected magnetic field. The diffusion region is localized so that 
the Petschek-type reconnection takes place. The slow-mode MHD shock is therefore formed, 
i.e., the gas pressure increases and the magnetic pressure decreases behind the shock front. The 
velocity of the reconnection outflow is the order of Alfven velocity. These are equivalent to the 
magnetic energy being converted to thermal and kinetic energy through the shock front. 

Figures 2c and 2d show the distribution in the case of 6' = 140° and 6 = 90°. The y- 
components of the initial magnetic fleld in these cases are less than that in the 6 = 180° case 
(the difference in the number of the contour indicates this situation) while the z-components 
are not equal to zero in these cases. Figure 3 shows the distribution of the x-component of 
the velocity normalized by waoCos(90° — 9/2) = vaq^cG on a line where y = at t/to = 20. The 
triangles display the results of the 9 = 180° case, the squares display the results of the 9 = 140° 
case, and the plus signs the 9 = 90° case. The inflow velocities Vin, especially near the current 
sheet (|a;|/Lo < 2), are proportional to t'Ao,cff- 

From the linearized MHD equations, the disturbance of velocity and magnetic fleld due 
to the Alfven wave is perpendicular to the initial (non-perturbed) magnetic fleld, Bq, and the 
wave vector, k (i.e., parallel to x Bq). In our case, the initial magnetic fleld is on the yz-p\ane. 
For simplicity, we assume that the wave vector is approximated to be parallel to the x (or —x) 
direction (the Alfven wave front is on the yz-plane). In this case, the perturbed components 
of velocity and magnetic fleld are on the yz-plane (no x-component) and perpendicular to 
the initial fleld. Hereafter, we call these velocity and magnetic fleld disturbances due to the 
Alfven wave as v± and B±. In this case, the perturbed components are described as v± = 
Vy sin (j) — V z cos (j) and Bj^ = By sin(f) — cos (p (see Figure 1). Note that the non-linear effect 
of shock propagating through the magnetic reconnection is included. However, it is difficult 
to remove the non-linear effect from these components. Therefore, this assumption would be 
acceptable for our flrst attempt to study the magnetic disturbance of non-perfectly antiparallel 
magnetic reconnection. 

In the case where there exists a magnetic fleld component perpendicular to the xy- 
plane, the fleld component perpendicular to the initial fleld, Bj_, is generated as a result of 
the reconnection although B± is not generated in the perfectly antiparallel reconnection (in the 
perfectly antiparallel reconnection case, nieans B^). The generated 5^ propagates as Alfven 
waves along fleld lines. The left four panels of Figures 4 and 5 show the time evolutions of —B± 
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in the cases of 6 = 140° and 6 = 90° respectively. The contour hnes show the magnetic field lines 
and the arrows show the velocity. The right panels of Figures 4 and 5 show the Bz distribution 
at t/to = 25. Following a certain magnetic field line, e.g., from the point {x/ L^.y / Lq) = (—3,20) 
in the 6 = 140° case, Bz once increases and then decreases with keeping its sign. This feature 
is consistent with the analytical solution in Petschek & Thorne (1967) (see Figure 3 in that 
paper). The point where Bz increases corresponds to Alfven wave and the point where Bz 
decreases corresponds to slow-mode MHD shock. 

3.2. Energy Fluxes of Alfven Waves and Magneto- acoustic Waves 

In this study, we focus on the disturbance energy fiuxes carried by the Alfven waves 
(-^Alfven) and magneto-acoustic waves (-Fsound) generated by the magnetic reconnection. These 
fluxes are measured on a line where y/Lo = ±10. 

i^Alfven = ±1-^ / -B±V±B\\y dx, (14) 

\ r ^'max 

Fsonnd = ±— / 6p V\\y dx, (15) 

where 5 means the difference from the initial value, Lx{= x^a.x — x^i^ = 26.2Lo) is the size of the 
computational domain in the x-direction. v\\y {Bi^\y) is the y-component of velocity (magnetic 
field) parallel to the initial magnetic field, Bq. v\\y and B^\y are described as 

VyByO^ + VzByoBzO , . 

ny= R 2 , R 2 ' (16) 

ByByO^ + BzBypBzO , , 

B\\y = R 2 , R 2 • (1") 

^yO + -D2O 

Here, ByQ (Bzq) is the ?/(2)-component of the initial magnetic field. The plus sign is on the 
y/Lo = 10 line and the minus sign on the y / Lq = —10 line. We also measure the time integral 
of each fiux: 

-E'Alfven = Lx / -^Alfven dt, (18) 







-E'Sound — Lx i -fsound dt. (19) 

-^Alfven and -Fsound are normalized by PoC'lo- -^Aifvcn and -Esound are normalized by PqC^qLI. 
The origin of the energies transported by the Alfven waves and magneto-acoustic waves is the 
magnetic energy released by the reconnection so that the magnetic energy in the simulation 
domain is also measured. 

Figures 6a-c show the time evolutions of the released magnetic energy (AE'mg = -Emgo ~ 
i^ing) in the simulation domain, where E^^ is defined by 

^mg = — / / 1^1' dx dy, (20) 

and -Emg is normalized by poC*so-^o- H^re, -Emgo is the initial magnetic energy in the simulation 
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domain. It is clear that the amount of the released magnetic energy decreases as the angle 
of magnetic shear, 9, decreases and the plasma /3 increases. In the final evolutional stage of 
each calculation, the magnetic energy in the simulation domain decreases linearly with time. 
Figure 6d shows the value of | dE^g/dt \ in the case of 6* = 180°, 140° and 90°. The magnetic 
energy released by the magnetic reconnection is equivalent to the Poynting flux entering into 
the reconnection region. Since the reconnection inflow, t^in, is expressed as v-m = €Va (x B in the 
case of Petschek-type reconnection, where e is the reconnection rate and roughly independent 
of B (e.g., Magara et al. 1996, Yokoyama & Shibata 1997), 

^^-^ 2^lv,^ocB'ocr''^\ (21) 



dt in 

where / is the size of reconnection region (Tanuma et al. 1999). The solid line in Figure 6d 
shows a line where | dE^g/dt |oc (3~^^^. It is clear from Figure 6d that the Equation (21) is 
almost satisfied independent of 9. 

Figures 7a and 7b show the time evolutions of the energy fluxes carried by the Alfven 
waves and magneto-acoustic waves, -FAifven and -Fsound, in the case of /3 = 0.1. If 6* = 180°, -FAifven 
is always zero because the perturbation perpendicular to the initial field is not generated. The 
time evolutions of the magneto- acoustic wave flux are similar in each case except when 9 < 60°. 
Those of the Alfven wave flux are also similar except for the cases where 9 = 180° and 9 < 60°. 
The magneto-acoustic and Alfven wave fluxes decrease after reaching the peak value. 

Figure 7c and 7d are the case of /3 = 1. The time evolutions of -FAifvcn and -Fsound in this 
case are similar in each 9 case respectively and similar to those in the /3 = 0.1 case. There is a 
following tendency; -FAifvcn at the late stage once increases and then decreases as 9 decreases, 
while -Fsound at the late stage decreases monotonically. Figures 7e and 7f show the results for 
the case of (3 = 20. Though the same features are seen, -FAifven and -Fsound are different by more 
than one order of magnitude when compared with the /3 = 1 case (the time at the late stage of 
the /3 = 20 case is about 3.5 times larger than that of the /3 = 1 case, i.e., the time integral is 
made about 3.5 times longer). 

4. Discussions and Summary 

Figure 8 shows the ratio (percentage) of the energies carried by the Alfven waves and 
magneto-acoustic waves (-^Alfven and -Esound) to the released magnetic energy at the final stage 
of simulation. In the case oi 9 = 180°, -^Alfven is equal to independent of the plasma (3 
because the magnetic field component perpendicular to the initial field is not generated by 
the reconnection. As 9 decreases the percentage of the energy carried by the Alfven waves 
increases up to 38.9%, 36.0%, and 29.5% in the case of /3 = 0.1, 1, and 20 respectively, and then 
decreases. The percentage is maximum when 9 = 80°. The percentage of the energy carried by 
the magneto- acoustic waves is roughly constant when 9 is relatively large in the cases of /3 = 0.1 
and 1. In the case of (3 = 20, it gradually decreases in decreasing 9. The maximum values are 

9 



16.2%, 25.9%, and 75.0% in the case of /3 = 0.1, 1, and 20 respectively. 

While the percentage of the energy carried by the Alfven waves is almost independent 
of /3, that of the energy carried by the magneto- acoustic waves changes by some factor. The 
ratio of the energy carried by the magneto-acoustic waves to the released magnetic energy 
becomes larger as /3 becomes larger in ^ > 40. This means that the significant part of the 
energy released by the magnetic reconnection is transported as a perturbation of gas pressure 
in a high-/9 plasma case. 

This result can give a suggestion to the study of high-/3 plasma astrophysical objects, 
e.g., the accretion disk. The accretion disk is thought to be weakly magnetized. The mag- 
netic reconnection is induced by the magnetorotational instability (e.g., Sano & Inutsuka 2001, 
Machida & Matsumoto 2003). Sano & Inutsuka (2001) showed that the heating rate is strongly 
related to the turbulent shear stress, which determines the efficiency of angular momentum 
transport. Therefore, the study of the magnetic reconnection in a high-/? plasma is important. 

The total non-radiative energy input to the solar coronal hole was estimated at 5 x 10^ 
ergs cm~^ s~^ by Withbroe (1988). For the acceleration of high-speed solar wind, some 1 x 10^ 
ergs cm~^ s~^ is required to be deposited at distances of several solar radii (see, e.g., Parker 
1991, and references therein). If the solar wind is accelerated by the energy flux of Alfven 
waves, this means that 20% of the energy released by reconnection events in the solar corona 
is transferred as a form of Alfven wave. Our results show that the energy larger than the 
required can be carried by the Alfven wave independent of /3 around the parametric region 
of 60° < 6 < 110°. These energies are converted to thermal energy through the dissipation of 
Alfven waves. 

The guide fleld, in this paper, can not contribute to the reconnection because of the 
2.5D approximation. Besides this, the reconnection progresses typically with the Alfven time 
scale, which depends on the magnetic fleld on the xy-plane (more exactly speaking, depend on 
the reconnection component of the Alfven velocity). It is therefore interesting that how the 
results change when the normalizations are changed: Time is re-normalized by the effective 
Alfven time (tAo,cff = Lo/vao,cs)- This means that the same time in Figure 6 is not the same 
time in Figure 9 because the effective Alfven velocity is different according to 6. The magnetic 
fleld is also re- normalized by the initial magnetic fleld which can reconnect (i?o sin 6*72). When 
^ = 0, the magnetic fleld is uniform and the reconnection does not take place. 

Figure 9 shows the difference of the re-normalized magnetic energy from the initial value 
as the function of re-normalized time in each 6 case. In the high-/? case, the magnetic energy is 
released at almost the same rate in the effective Alfven time when 6 > 40°. The lower /3 is, the 
larger the discrepancy of the release rate becomes. Figure 10 shows the re-normalized energy 
fluxes and the time integral of those as the function of re-normalized time. The energy flux and 
its time integral of Alfven waves show the time variation similar to each other independent of 
p. Those of magneto-acoustic waves also show the time variation similar to each other except 
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when 9 is relatively small, although there is a bit of j3 dependence. Compared with Figure 7, 
the peak positions are matched. 

Figure 11 indicates the percentage of the energies carried by the Alfven waves and 
magneto-acoustic waves to the released magnetic energy. The vertical axis is re-normalized as 
same as the above-mentioned way. This shows clearly that the amount of the energy carried 
by the Alfven waves has almost the same dependence on 9 independent of /3 if the energy and 
time are scaled by the effective magnetic energy and Alfven time. This is equivalent to that the 
magnetic configuration is important rather than the field strength relative to the gas pressure 
for the energy release rate in the effective Alfven time. On the other hand, the amount of the 
energy carried by the magneto-acoustic waves shows the (3 dependence. 

In this paper, we have reported the results for 2.5D MHD simulations of the magnetic 
reconnection. When magnetic fields are non-perfectly antiparallel, the magnetic field compo- 
nent perpendicular to the initial field is generated and propagates as the Alfven wave. We have 
measured the energy fluxes of Alfven waves and magneto-acoustic waves. Magneto-acoustic 
waves are related to fast mode waves in the high-/? case {(3 = 20) and slow mode waves in 
the low-/9 case (/3 = 0.1). The energy carried by the Alfven waves is more than 30% of the 
energy released by the magnetic reconnection at the maximum. That value satisfies the re- 
quirement for energy flux of Alfven waves necessary for acceleration of high-speed solar wind 
in the nanofiare coronal heating model. For more exact discussion, 3-dimensional simulations 
with more realistic plasma distribution are necessary. 
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Table 1. Values of Parameters in Simulations 



plasma /3 


Alfven velocity, vao 




0.1 


-3.5 


- 2500 
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- 1.1 


-800 


20 


-0.24 


- 180 
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(□) t/t(j=0.0 (b) t/t(j=20.0 (c) t/to=20.0 (d) t/t(j=2D.O 

e=iso° e=iso' 6=^40' 8= 90° 




Fig. 2. Distribution of the logarithmic pressure in the case of (a) 6 — 180° at initial state, (b) 6 = 180° 
at t/to = 20.0, (c) = 140° at t/to = 20.0 and (d) ^ 90° at t/to = 20.0. The plasma (3 is 0.1. The black 
lines show the magnetic field lines, and the arrows show the velocity. 
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Fig. 3. The distribution of the x-component of the velocity normalized by uaoCos(90° — 9/2) = i^Ao.cff on 
a line where y = at t/to = 20. The triangles display the results of the 9 = 180° case, the squares display 
the results of the 9 — 140° case, and the plus signs display the results of the 9 = 90° case. The plasma 
/? is 0.1. The inflow velocities, in this case, are roughly proportional to cos (90° — 9/2). This can be 
interpreted as Vx oc WAo.off , where UAo.off is the Alfvcn velocity in the computational plane. 
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Fig. 4. The left four panels show the time evolution of the magnetic field component perpendicular to 
the initial field, ~B±, in the case oi 9 = 140°. The right panel shows the distribution at t/to = 25. The 
plasma /3 is 0.1. The black lines show the magnetic field lines, and the arrows show the velocity. 
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Fig. 5. The left four panels show the time evolution of the magnetic field component perpendicular to 
the initial field, —B±, in the case of 6* = 90°. The right panel shows the B, distribution at t/to = 25. The 
plasma /3 is 0.1. The black lines show the magnetic field lines, and the arrows show the velocity. 
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Fig. 6. (a)-(c) Time evolution of the released magnetic energy in the computational domain in each 9 
case. The plasma /? are (a) 0.1, (b) 1, and (c) 20. (d) The dependence of the magnetic energy release 
rate in the late stage on /3 in each 9 case. The solid line shows a line where | dE^^g/dt |oc The 
magnetic energy release rate is proportional to f3~^^^, which is consistent with the theoretical prediction, 
independent of 9. 
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Fig. 7. Time evolution of tlie energy fluxes, F ^ (a, c, e) Alfven waves and (b, d, f) magneto-acoustic 
waves in each Q case. The plasma /3 are (a, b) 0.1, (c, d) 1.0 and (e, f) 20.0. 
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Fig. 8. The percentage of the energies carried by the Alfven waves (i?Aifvcn/ Aiding) and magneto-acoustic 
waves (-Esound/Ai?mg) to the released magnetic energy at the final stage of simulation in the case of 
61 = 180°. 
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Fig. 9. Time evolution of the magnetic energy in the computational domain in each 9 case. Compared 
with Figure 6, the horizontal axis, time, is re- normalized by the effective Alfven time (tAo,off = ^o/'^Ao,cff)- 
The vertical axis, the difference of the magnetic energy from the initial value, is also re-normalized as 
-A£;* = -I\E,-^^/{BQ^sin^e/2). The line styles are the same as Figure 6. 
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Fig. 10. Time evolution of tlie re-normalized energy fluxes, F*, of (a, c, e) Alfven waves and (b, d, f) 
magneto-acoustic waves in each Q case. The energy fluxes are re-normalized as F* = F j (uAo,cff-Bo^sin^0/2). 
The plasma /3 are (a, b) 0.1, (c, d) 1.0 and (e, f) 20.0. The line styles are the same as Figure 7. 
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Fig. 11. The percentage of the energies carried by tlie Alfven waves (iJ^jf^^j^/AiJjJjg) and magneto-acous- 
tic waves {Ei^oundl ^-^i-ag) to the released magnetic energy. The released magnetic energy is re-normalized 
by the energy of the initial magnetic field which can reconnect. The symbols are the same as Figure 8. 
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